Quantum information processing with circuit quantum electrodynamics 
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We theoretically study single and two-qubit dynamics in the circuit QED architecture. We focus 
on the current experimental design [Wallraff et al, Nature 431, 162 (2004); Schuster et at, Nature 
445, 515 (2007)] in which superconducting charge qubits are capacitively coupled to a single high- 
Q superconducting coplanar resonator. In this system, logical gates are realized by driving the 
resonator with microwave fields. Advantages of this architecture are that it allows for multi-qubit 
gates between non-nearest qubits and for the realization of gates in parallel, opening the possibility 
of fault-tolerant quantum computation with superconduting circuits. In this paper, we focus on one 
and two-qubit gates that do not require moving away from the charge-degeneracy 'sweet spot'. This 
is advantageous as it helps to increase the qubit dephasing time and does not require modification 
of the original circuit QED. However these gates can, in some cases, be slower than those that do 
not use this constraint. Five types of two-qubit gates are discussed, these include gates based on 
virtual photons, real excitation of the resonator and a gate based on the geometric phase. We also 
point out the importance of selection rules when working at the charge degeneracy point. 

PACS numbers: 03.67.Lx, 73.23. Hk, 74.50.+r, 32.80.-t 
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I. INTRODUCTION 

Superconducting circuits based on Josephson junc- 
tions P, are currently the most experimentally ad- 
vanced solid-state qubits. The quantum behavior of 
these circuits has been experimentally tested at the 
level of a single qubit 0, 0, @, @, 0] and of a pair of 
qubits [1, H, E3, EH, EH- The first quantitative experi- 
mental study of entanglement in a pair of coupled super- 
conducting qubits was recently reported [l3[ . 

In this paper, we theoretically study quantum infor- 
mation processing for superconducting charge qubits in 
circuit QED % El, El EE EE El , focusing on two-qubit 
gates. In this system qubits are coupled to a high Q 
transmission line resonator which acts as a quantum bus. 
Coupling of superconducting qubits through a quantum 
bus has already been studied by several authors and in 
different settings. In particular, cou plin g using a lumped 



21 , |2l |2JjM gajman extended 
23, EE El SHil, a current- 



LC oscillator [19|, [20j 
ID or 3D resonator 



biased Josephs on j unction acting as an anharmonic os- 
cillator 32 , 34 , [H, [H, H3] or using a mechanical oscilla- 



tor |3J, 13J, |40| were studied. Here, we focus on circuit 
QED with charge qubits [l4[ and consider the constraints 
of the current experimental design [!, El, El, E3 ■ As wc 
will show, while this architecture is simple, it allows for 
many different types of qubit-qubit interactions. These 
gates have the advantage that they can be realized be- 
tween non-nearest qubits, possibly spatially separated by 
several millimeters. In addition to being interesting from 
a fundamental point of view, this is highly advantageous 
in reducing the complexity of multi-qubit algorithms [4l| . 
Moreover, it also helps in reducing the error threshold 



required for reaching fault-tolerant quantum computa- 
tion (42j. Furthermore, some of the gates that will be 
presented allow for parallel operations (i.e. multiple one 
and two-qubit gates acting simultaneously on different 
pairs of qubits). This feature is in fact a requirement 
for a fault-tolerance threshold to exist [43[ , and this puts 
circuit QED on the path for scalable quantum computa- 
tion. 

Another aspect addressed in this paper is the 'quality' 
of realistic implementations of these gates. To quantify 
this quality, several measures, like the fidelity, have been 
proposed j44[ . A fair evaluation and comparison of these 
measures for the different gates however requires exten- 
sive numerical calculations including realistic sources of 
imperfections and optimization of the gate parameters. 
In this work, we will rather present estimates for the 
quality factor [f| of the gates as obtained from analytical 
calculations. Initial numerical calculations have showed 
that, in most situations, better results than predicted by 
the analytical estimates can be obtained by optimization. 
The quality factors presented here should thus be viewed 
as lower bounds on what can be achieved in practice. 

Five types of two-qubit gates will be presented. First, 
we discuss in section IIVI gates that are based on tuning 
the transition frequency of the qubits in and out of reso- 
nance with the resonator by using dc charge or flux bias. 
As will be discussed, this approach is advantageous be- 
cause it yields the fastest gates, whose rate given by the 
resonator-qubit coupling frequency. A problem with this 
simple approach is that it takes the qubits out of their 
charge-degeneracy 'sweet spot', which can lead to a sub- 
stantial increase of their dephasing rates ||. Moreover, 
changing the qubit transition frequency over a wide range 
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of frequencies can be problematic if the frequency sweep 
crosses environmental resonances j4^| . 

To address these problems, we will focus in this pa- 
per on gates that do not require dc excursions away from 
the sweet spot. Requiring that there is no dc bias is a 
stringent constraint and the gates that are obtained will 
typically be slower. However, the resulting gate qual- 
ity factor can be larger because of the important gain 
in dephasing time. The first of these types of gates rely 
on virtual excitation of the resonator (section PVTl . This 
type of approach was also discussed in Refs. [lj, [29| and 
here we will present various mechanisms to tune this type 
of interaction. The dispersive regime that is the basis of 
the schemes relying on virtual excitations of the resonator 
can also be used to create probabilistic entanglement due 
to measurement. This is discussed in section fVTl Next, 
we consider gates that are based on real photon popu- 
lation of the resonator (section IVII[) . For these gates, 
selection rules will set some constraints on the transi- 
tions that can be used. Finally, we discuss a gate based 
on the geometric phase which was first introduced in the 
context of ion trap quantum computation prl l47j . 

Before moving to two-qubit gates, we begin in sec- 
tion [IT] with a brief review of circuit QED and, in sec- 
tion |TTT1 with a discussion of singlc-qubit gates. A table 
summarizing the expected rates and quality factors for 
the different gates is presented in the concluding section. 



II. CIRCUIT QED 

A. Jaynes-Cummings interaction 

In this section, we briefly review the circuit QED ar- 
chitecture first introduced in Ref. [lij and experimentally 
studied in Refs. 0, [IH, [l7| ■ Measurement-induced de- 
phasing was theoretically studied in Ref. [l8[. As shown 
in Fig.rn this system consists of a superconducting charge 



qubit 



strongly coupled to a transmission line 
Near its resonance frequency u r , the 



resonator [5( 

transmission line resonator can be modeled as a simple 
harmonic oscillator composed of the parallel combination 
of an inductor L and a capacitor C . Introducing the an- 
nihilation (creation) operator aS^\ the resonator can then 
be described by the Hamiltonian 



uj r a'a, 



(2.1) 



with u r = 1/vLC and where we have taken % = 1. Using 
this simple model, the voltage across the LC circuit (or, 
cquivalently, on the center conductor of the resonator) 
can be written as Vlc = Kms( at + a )> where V^ ns = 
^Jhui r /2C is the rms value of the voltage in the ground 
state. An important advantage of this architecture is the 
extremely small separation 6^5 /im between the center 
conductor of the resonator and its ground planes. This 
leads to a large rms value of the electric field E® = 
v r°ms/ b ~ °- 2 v / m for typical realizations 0, QJ, M, 03 ■ 




FIG. 1: (Color online) Layout and lumped element version of 
circuit QED. Two superconducting charge qubit (green) are 
fabricated inside the superconducting ID transmission line 
resonator (blue). 



Multiple superconducting charge qubits can be fabri- 
cated in the space between the center conductor and the 
ground planes of the resonator. As shown in Fig. [U we 
will consider the case of two qubits fabricated at the two 
ends of the resonator. These qubits are sufficiently far 
apart that the direct qubit-qubit capacitance is negligi- 
ble. Direct capacitive coupling of qubits fabricated inside 
a resonator was discussed in Ref. [29[. An advantage of 
placing the qubits at the ends of the resonator is the 
finite capacitive coupling between each qubit and the in- 
put or output port of the resonator. This can be used to 
independently dc-bias the qubits at their charge degen- 
eracy point. The size of the direct capacitance must be 
chosen in such a way as to limit energy relaxation and 
dephasing due to noise at the input/output ports. Some 
of the noise is however still filtered by the high-Q res- 
onator [l4j]. We note, that recent design advances have 
also raised the possibility of eliminating the need for dc 
bias altogether. [I?) 

In the two-state approximation, the Hamiltonian of the 
jth qubit takes the form 



H n = — 



(2.2) 



where E e \. = 4£'c j (l — 2n Sj ) is the electrostatic energy 
and -Ejj = £;™ ax cos(7r$j/$o) is the Josephson coupling 
energy. Here, Ec, = e 2 /2Cs j is the charging energy with 
Csj the total box capacitance. n gj = C gj V gj /2e is the 
dimensionless gate charge with C gj the gate capacitance 
and V gj the gate voltage. iJ™ ax is the maximum Joseph- 
son energy and $j the externally applied flux with $o 
the flux quantum. Throughout this paper, the j sub- 
script will be used to distinguish the different qubits and 
their parameters. 

With both qubits fabricated close to the ends of the 
resonator (antinodes of the voltage), the coupling to the 
resonator is maximized for both qubits. This coupling 
is capacitive and determined by the gate voltage V gj = 
V g c + Vlc , which contains both the dc contribution V^ c 
(coming from a dc bias applied to the input port of the 
resonator) and a quantum part Vlc- Following Ref. [lij], 
the Hamiltonian of the circuit of Fig. Q] in the basis of 
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the cigenstates of H q . takes the form 

TT t ^7 

H =u) r a>a+ ~ az i 

3 = 1,2 

~ E Sj - c j°z 3 + Sja Xj ) (at + a) 

.7=1,2 



(2.3) 



where u> aj = yj E\ + [4Ec 3 (1 — 2n g _j)] 2 is the transi- 
tion frequency of qubit j and gj = e(C g j /Csj)V r ° ms /h 
is the coupling strength of the resonator to qubit j. 
For simplicity of notation, we have also defined pj 



1 - 2n 



9,31 



COS( 



and s, = sin i 



where 9j 



arctan(£'j j /Eg, (1 — 2n g j)) is the mixing angle [14(. 

When working at the charge degeneracy point = 
1/2, where dephasing is minimized [5j , and neglecting fast 
oscillating terms using the rotating-wave approximation 
(RWA) , the above resonator plus qubit Hamiltonian takes 
the usual Jaynes-Cummings form fn\ 



E 

3=1,2 



E 9 3 

3=1,2 



(2.4) 

This coherent coupling between a single qubit and the 
resonator was investigated experimentally in Ref. [H, [TEl . 
[H,[l3|- In particular, in Ref. [|| high fidelity single qubit 
rotations were demonstrated. 



B. Damping 

Coupling to additional uncontrollable degrees of free- 
dom leads to energy relaxation and dephasing in the sys- 
tem. In the Born-Markov approximation, this can be 
characterized by a photon leakage rate k for the res- 
onator, an energy relaxation rate 71 t j and a pure de- 
phasing rate 7^ for each qubit. In the presence of these 
processes, the state of the qubit plus cavity system is de- 
scribed by a mixed state p(t) whose evolution follows the 
master equation [5l| 



p = -i[H,p] + nV{a]p + H,3 V \°- ] ]P 

3=1,2 

+ E y^KlA 

J=l,2 



(2.5) 



where V[L]p = {2Lptf - L^Lp - pL^h) /2 describes the 
effect of the baths on the system. 



C. Typical system parameters 

In this section, we give realistic system parameters. 
The resonator frequency Lo r /2iT will be assumed to be 
between 5 and 10 GHz. The qubit transition frequencies 
can be chosen anywhere between about 5 to 15 GHz, and 
are tunable by applying a flux though the qubit loop. In 



the schematic circuit of Fig. [T] both qubits arc affected by 
the externally applied field, but the effect on each qubit 
will depend on the qubit's loop area. Coupling strengths 
g/2ir between 5.8 and 100 MHz have been realized exper- 
imentally [HI, [ljj and couplings up to 200 MHz should 
be feasible. 

Rabi frequencies of 50 MHz where obtained with a 
sample of moderate coupling strength g/2n = 17 MHz Q 
and an improvement by at least a factor of two is realistic. 

The cavity damping rate k is chosen at fabrication time 
by tuning the coupling capacitance between the resonator 
center line and it's input and output ports. Quality fac- 
tors up to Q ~ 10 6 have been reported for under-coupled 
resonators [5(| , corresponding to a low damping rate 
k/2tt = <jj r /2nQ ~ 5 KHz for a uj r /2ir = 5 GHz resonator. 
This results in a long photon lifetime 1/k of 31 ps. To 
allow for fast measurement, the coupled quality factor 
can also be reduced by two or more orders of magnitude. 

Relaxation and dephasing of a qubit in one realization 
of this system were measured in Ref. @. There, T\ = 
7.3 ps and T2 = 500 ns were reported. These translate to 
7i/2tt = 0.02 MHz and 7 /2tt = ( 72 - 7i/2)/2tt = 0.31 
MHz. 



III. 1-QUBIT GATES 

Single qubit gates are realized by pulses of microwaves 
on the input port of the resonator. Depending on the 
frequency, phase and amplitude of the drive, different 
logical operations can be realized. External driving of 
the resonator can be described by the Hamiltonian 



e* k {t)ae 



(3.1) 



where s k (t) is the amplitude and LOd k the frequency of 
the kth external drive. Throughout this paper, the k 
subscript will be used to distinguish between the different 
drives and the drive-dependent parameters. 

For simplicity of notation, we first consider the situa- 
tion where there is a single qubit and drive present. We 
will also assume that the qubit is biased at its optimal 
point and use the RWA. The Hamiltonian describing this 
situation is H = H jc + Hd with j = k = 1 . 

Logical gates are realized with microwaves pulses that 
are substantially detuned from the resonator frequency. 
With a high-Q resonator, this means that a large frac- 
tion of the photons will be reflected at the input port. To 
get useful gate rates, we thus work with large amplitude 
driving fields. In this situation, quantum fluctuations in 
the drive are very small with respect to the drive am- 
plitude and the drive can be considered, for all practical 
purposes, as a classical field. In this case, it is convenient 
to displace the field operators using the time-dependent 
displacement operator [53[ 



D(a) = cxp 



aa' 



a a) 



(3-2) 
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Under this transformation, the field a goes toa+a where 
a is a c-number representing the classical part of the field. 
The displaced Hamiltonian reads 

H = D f (a)HD(a) - iD\a)D{a) 

= uj r a) a + ~^ a z ~ 9 {p^cr- + u+a) (3-3) 
-g(a*a- + aa+) , 

where we have chosen a(t) to satisfy 

a = -i(j r a - ie(t)e- iUdt . (3.4) 

This choice of a is made so as to eliminate the direct drive 
on the resonator Eq. (|3.1[) from the effective Hamiltonian. 

In the case where the drive amplitude e is independent 
of time, and by moving to a frame rotating at the fre- 
quency LDd for both the qubit and the field operators, we 
get 

H = A r a^a H — 7f~ a z ~ 9 ( a ^ (J - + <J + a ) H — ^~° x ' (^*^) 

where we have dropped any transient in a(t). In the 
above expression, we have defined A r = uj r — LUd which is 
the detuning of the cavity from the drive, A a = uj a — cu c [ 
the detuning of the qubit transition frequency from the 
drive and flu is the Rabi frequency: 



(3.6) 



In the limit where A,, is large compared with the res- 
onator half-width re/2, the average photon number in the 
resonator can be written as n w (e/A r ) 2 . In this case, 
the Rabi frequency takes the simple form flu m 2g\ffi 
expected from the Jaynes-Cummings model. 

We note that the effect of damping can be taken into 
account by performing the transformation (13.21) on the 
master equation (|2.5|) rather than on Schrodinger's equa- 
tion. For completeness, this is done in appendix [XI Since 
in this paper we are interested in the qubit dynamics un- 
der coherent control rather than measurement, we will 
be working in the regime where A r > re and as such 
can safely ignore the effect of re on f2#. For a detailed 
discussion of measurement in this system, see Ref. [l8j . 



A. On-resonance: bit-flip gate 

For quantum information processing, it is more ad- 
vantageous to work in the dispersive regime where A = 
uj a — Lu r ) is much bigger than the coupling g. One advan- 
tage of this regime is that the pulses aimed at controlling 
the qubit are far detuned from the resonator frequency 
and are thus not limited in speed by its high quality fac- 
tor. Another advantage is that the high quality resonator 
filters noise at the far detuned qubit transition frequency 
and effectively enhances the qubit lifetime (l4| . 



To take into account that we are working in this dis- 
persive regime, we eliminate the direct qubit-resonator 
coupling by using the transformation 



U = exp 



9 , + 

— (o'er- - acr+, 



(3.7) 



Using the Hausdorff expansion to second order in the 
small parameter A = g/A 



A 



e - AA He AX = H + X[H,X] + — [[J^ X\ , X 

with X = (a'cr_ — acr + ), yields [3 



(3.8) 



H x w A r a)a + - [A a + 2 X 

A„ rt R 



a'a+- 



A t i a i 



(3.9) 



where we have defined \ = 5 2 /A and A a = Q a — ujd 
with u> a — uj a + x- Since the resonator is driven far 
from the frequency band uj r ± \ where cavity population 
can be large, we have that (a'a) ~ (this is because we 
are working in a displaced frame with respect to the res- 
onator field). As a result, we have therefore dropped the 
ac-Stark shift in the second line of the above expression. 

By choosing A a = 0, the above Hamiltonian generates 
rotations around the x axis at a rate These Rabi 

oscillations have already been observed experimentally 
in circuit QED with close to unit visibility [3j . Changing 
A a , flu and the phase of the drive can be used to rotate 
the qubit around any axis on the Bloch sphere [54| . 

In the situation where many qubits are fabricated in 
the resonator and have different transition frequencies, 
the qubits can be individually addressed by tuning the 
frequency of the drive accordingly. It should therefore 
be possible to individually control several qubits in the 
circuit QED architecture. 



B. Off-resonance: phase gate 

It is useful to consider the situation where the drive is 
sufficiently detuned from the qubit that it cannot induce 
transitions, but is of large enough amplitude to signifi- 
cantly ac-Stark shift the qubit transition frequency due 
to virtual transitions. To obtain an effective Hamiltonian 
describing this situation, we start by adiabatically elim- 
inating the effect of direct transitions of the qubit due 
to the drive. This is be done by using on Eq. (|3.5[) the 
transformation 



U = exp[f3*a+ - /3a-] 



(3.10) 



to second order in the small parameter (3 = A a . In 
a second step, we again take into account the fact that 
the qubit is only dispersively coupled to the resonator 
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by using the transformation of Eq. (|3 . T[) to second order. 
These two sequential transformation yield 



At 1 fx I^r 



(3.11) 



The last term in the parenthesis is an off-resonant ac- 
Stark shift caused by virtual transitions of the qubit. 
This shift can be used to realize controlled rotation of 
the qubit about the z axis. The rate of this gate can 
be written in terms of the average photon number n in- 
side the resonator as ~ 2gV^ X (f2fl/2A a ). To get fast 
rotations, one must therefore choose large values of the 
coupling constant g and large n while keeping the ratio 
Qfi/2A a small to prevent real transitions. 

Finally, it is important to point out that in the sit- 
uation where multiple qubits are present inside the res- 
onator, each qubit will suffer a frequency shift when other 
qubits are driven. These frequency shifts will have to be 
taken into account or canceled by additional drives. 



Coherent control vs. measurement 



-2jt 10 MHz 




FIG. 2: (Color online) Resonator transmission (dashed lines) 
and corresponding phase shifts (full lines) for the two qubit 
states (blue: ground, red: excited). The numbers are cal- 
culated using g/2n ~ 100 MHz and g/A ~ 0.1. At the far 
detuned frequencies required for single qubit gates, the phase 
shift of the resonator field does not depend on the qubit state. 
In this case, there is negligible entanglement between the res- 
onator and the qubit. 

As mentioned above, in the dispersive regime, driving 
the cavity close to its resonance frequency leads to a mea- 
surement of the qubits. As discussed in Refs. [3, HI], this 
is due to entanglement of the qubit with the resonator 
field generated by the term x a<[aa z of Eq. (|3.9|) . Indeed, 
because of this term, the resonator frequency is pulled to 
u> r ± x depending on the state of the qubit. The possi- 
ble resonator transmissions, corresponding to the qubit 
in the ground (blue) or excited (red) state, are shown 
(dashed lines) along with the corresponding phase shifts 
(full lines) in Fig. d 

As is seen from this figure, only around uj r ± x is 
there significant phase shift and/or resonator transmis- 
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FIG. 3: (Color online) Quality factor of the phase gate, with 
respect to the measurement-induced dephasing rate F m (blue, 
full line) and with respect to 72 +F m (red, dashed line) where 
I/72 = 500 ns. This is plotted vs. the detuning between 
the control drive and the qubit transition frequency. The 
amplitude e of the drive is chosen such that Q.R,/2A a =0.1 for 
all values of the detuning. The other parameters are g/2n = 
100 MHz, g/A = 0.1 and k/2tt = 0.5 MHz. The discontinuity 
in r m around (uj a —ui) /2n = A = 1000 MHz is due to the fact 
at that frequency, the drive is on resonance with the resonator 
and therefore corresponds to a measurement. This is not a 
region where a phase gate would be operated. 



sion change for the information rate about the qubit 's 
state to be large at the resonator output [13] . In other 
words, only around these frequencies is entanglement be- 
tween the resonator and the qubit significant. However, 
when coherently controlling the qubit using the flip and 
phase gates discussed above, the resonator is irradiated 
far from w r ± x- As shown in Fig. [21 since we are work- 
ing in the dispersive regime where |A| 3> g, there is no 
significant phase difference in the resonator output be- 
tween the two states of the qubit at these very detuned 
frequencies. As a result, there is no significant unwanted 
entanglement with the resonator when coherently con- 
trolling the qubit. 

An additional benefit of working at these largely de- 
tuned irradiation frequencies is that the resonator is only 
virtually populated and the speed of the gates is not lim- 
ited by the high Q of the resonator. These two aspects 
lead to high quality single qubit gates Q ■ 

The above discussion can be made more quantitative 
by introducing the rate T m of dephasing induced by the 
control drive (corresponding to measurement-induced de- 
phasing) [rsj : 



where 



n± 



Kx 2 (n + + n_) 
( K /2) 2 + A2+ x 2 



( K /2)2 + (A r ± X ) 2 



(3.12) 



(3.13) 



is the steady-state average photon number inside the res- 
onator for a qubit in the ground (— ) or excited (+) state. 
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In practice, this rate will always be much smaller than 
the intrinsic dephasing rate I/T2 of the qubit. For ex- 
ample, for the bit-flip gate, a Rabi rate of £Ir/2tt = 1 00 
MHz with g/2n = 100 MHz and g/A = 0.1 yields a 
measurement-induced dephasing time 1 /T m of the order 
of a few milliseconds. Clearly, this is not a limitation in 
practice. This is illustrated for the phase gate in Fig. [3] 
where the quality factor 



Q<fi 



n 2 R /2A a 

2T A 



(3.14) 



is plotted as a function of the detuning of the drive with 
respect to the qubit transition frequency. In this figure, 
the full blue line is the quality factor (Q m ) calculated 
using the measurement-induced dephasing rate r rra and 
the dashed red line the quality factor (Qt) using the total 
rate = 72 + T m assuming a dephasing time I/72 of 
500 ns. 

For the phase gate, a dephasing time of I/72 = 500 
ns with a rate of il^/AnA a = 40 MHz at a detuning 
(co a —Ud)/2TT = 2000 MHz yields a quality factor of ~ 60. 
For the bit-flip gate, a Rabi rate of 100 MHz yields a 
quality factor of ~ 157. 



D. AC-dither: phase gate 

Another approach to produce a single-qubit phase 
gate, is to take advantage of the quadratic dependence 
of the qubit transition frequency on the gate voltage (or 
flux) to shift the qubit transition frequency. This can 
be done by modulation of these control parameters at 
a frequency that is adiabatic with respect to the qubit 
transition frequency. 

Focussing on the single qubit Hamiltonian (12.21) . we 
i d 

a modulation of the gate voltage that is slow compared 
to the qubit transition frequency. In this situation, it 
is useful to move to the adiabatic basis. The relation 
between the original (<7j) and adiabatic (dj) basis Pauli 
operators is given by 




CO 

- 4— I 

'cz 

=> 

< 

CD 
CZ 
LU 



0.25 0.5 0.75 1 
gate voltage, n g 

FIG. 4: (Color online) Refocussing due to ac-dither. The 
qubit energy splitting is shown as a function of gate charge. 
Ac-dither acts as refocussing by sampling the positive and 
negative dependence of w„ (t) with n g . 



which should be compared to the bare value Ej. Here, 
we have dropped terms rotating at 2w ac and higher order 
in the dither amplitude. Voltage ac-dither can therefore 
be used to blue-shift the qubit transition frequency (flux 
dither around the flux sweet would cause a red shift). As 
will be discussed below, this can also be used to couple 
qubits when the dither frequency is larger than the cou- 
pling strength (but still slow with respect to the qubit 
transition frequency) . 

Because ac-dither acts as a continuous spin-echo [55j . 
it can be realized with minimal dephasing of the qubit. It 
therefore appears to be more advantageous than dc bias 
of the control parameters. For example, if the qubit is 



take n g (t) = n dc + n d (t), where n d (t) = nf sin(w ac t) is dc-biased An g away from the charge degeneracy point, 

n dc = 1/2 + An g , noise Sn g (t) in the bias will cause de- 
phasing due to fluctuations (duj a / ' dn g )5n g (t) in the qubit 
transition frequency. As illustrated in Fig. 0] if the dc- 
offset An g is small compared to the dither amplitude n^ c , 
then under dither, both signs of duj a /dn g will be probed 
leading to (partial) cancelation of the unwanted fluctu- 
ating phase. 

This cancelation can be seen more explicitly by assum- 
ing small excursions away from the charge degeneracy 



a z = cos <d(t)a z — sin <d(t)a x , 
a x = sin0(i)<72 + cosQ(t)a x , 



(3.15) 



where 9(t) = arctan[-E c i(0/-^'j]- In this basis, the qubit 
Hamiltonian reads 



point such that <d(t) w — A(An fl 



: sin o> ac i) , where we 



H = - 



(3.16) 



with 



have defined A = 8Eq/Ej. Assuming that the qubit is 
not too deep in the charge regime such that An^ c is small, 
we expand to first order in An g and n^ c to obtain 



{4£ c [l - 2n dc - 2n d (t)}} 2 the in- 
stantaneous splitting. Because of the quadratic depen- 
dence with gate charge, the average part of the qubit 
splitting is larger than its bare value. For example, set- 
ting n g = 1/2 and assuming small dither amplitude, we 
obtain 



H 



— o,. 



U) 



ad 



ifi]f K c ) : 



(3.17) 



+ AE c 5n g (t) [—\An g Jo (\n^) - 2 Ji(An^ c ) sinw ac t] ~a z 

+ 4E c 6n g (t) [J (Xn a g c ) - 2XAn g Ji(An^ c ) sin w ac i] 5 X , 

(3.18) 

where the J n (z) arc Bessel functions of the first kind. 
In this expression, we have dropped higher order Bessel 
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functions and have added the gate charge noise Sn g (t) to 
,(*)• 



The second term in Eq. (|3. 18|) . proportional to a z , 
leads to pure dephasing while the last term leads to 
mixing of the qubit. Focusing on pure dephasing, we 
obtain from the Golden rule 



1 1 /64££ 



Tj> 2 \ Ej 

"9 



&n 2 q S Sn (0) 



K c ) 2 



\ 2 

; ){ 



(3.19) 



where Ss ng (^ac) is the spectral density of the charge 
noise. 

For both the dc and ac bias, the dephasing rate 1/T^ 
increases linearly with the amplitude of the shift of the 
transition. However, blue shifting of the qubit using ac- 
dither n^ c produces less dephasing than using the static 
bias An g provided that the dither frequency is much 
higher than the characteristic frequency of the noise so 
that 



S Sn (±w ac ) < S Sng (0). 



(3.20) 



This approach should therefore efficiently protect the 
qubit from low frequency (i.e. 1//) noise. As- 
suming An g = and using Eq. (|3.17p . the qual- 
ity factor of this gate can be estimated as Q ~ 
(Ej/64El) 2 /{S Sng (-oj ac ) + S 5 „ 3 (+o> ac )). This type of 
stabilization of logical gate by ac-fields was also studied 
in Ref. (H. 



IV. DIRECT COUPLING BY VARIABLE 
DETUNING 

In the layout of Fig. [TJ qubit-qubit interaction must 
be mediated by the resonator. It is therefore reason- 
able to expect that the limiting rate on two-qubit gates 
be the coupling strength g. Gates at this rate can be 
implemented by taking the qubits in and out of reso- 
nance with the resonator frequency. When a qubit is 
far off resonance, it only dispcrsivcly couples to the res- 
onator through the coupling x°^ a(J z, where x = 9 2 /A. 
This interaction can be made small by working at large 
dctunings A. In this situation there is no significant 
qubit-resonator interaction. The interaction is turned 
on by tuning the qubit transition frequency back in res- 
onance with the resonator. In this case, vacuum Rabi 
flopping at the frequency 2g will entangle the qubit and 
the resonator. It is know from ion-trap quantum com- 
puting [57ll58|. and is further discussed in Section fVII D[ 
how to use this type of interaction to mediate qubit-qubit 
entanglement. 

Tuning of the qubit transition frequency could be real- 
ized by applying flux pulses through the individual qubit 
loop. This would require adding flux lines in proximity 
to the qubits. Voltage bias using individual bias lines 
could also be used, but this likely introduce more noise 



than flux bias. Moreover, in both cases, this will take 
the qubits away from their sweet spot, possibly increas- 
ing their dephasing rates Alternatively, Wallquist et 
al. have suggested that a similar tuning of A can be re- 
alized by fabricating a resonator whose frequency is itself 
tunable [27j |. While this is a promising idea, one draw- 
back is that any noise in the parameter controlling the 
resonator frequency will lead to dephasing of photon su- 
perpositions, lowering the expected gate quality. 

In addition to dc-bias, it is possible to use any of the rf 
approaches discussed in the previous section to tune the 
qubit transition frequency. Moreover, the FLICFORQ 
protocol [59| . discussed in the next section and in ap- 
pendix [C] could also be used. As shown in appendix O 
this would yield qubit-resonator coupling at the rate g/2. 
However, for these approaches to be useful here, very 
large rf amplitudes would be required to cover the large 
range of frequency needed to turn on and off the qubit- 
resonator interaction. This is especially true in the pres- 
ence of many qubits fabricated in the same resonator. 
A FLICFORQ-typc protocol [5^] was also suggested for 
flux qubits coupled to a LC oscillator in Ref. [lji[. In 
this case, the authors considered quantum computation 
in the basis of the qubit dressed by a rf-drive directly 
applied directly to the flux qubit. 

In summary, this type of gate relying on tuning of the 
qubit or resonator frequency is advantageous because it 
operates at the optimal rate g. However, it requires ei- 
ther additional bias lines and extra dephasing or large 
amplitude rf-pulses. In the next sections, we will focus 
on gates that do not require additional tuning but only 
rely on rf-drive of the resonator of more moderate am- 
plitudes. While these gates will be typically slower than 
the gates discussed here, they can be implemented with- 
out modification of the original circuit QED design and 
do not suffer from the above problems. Gates relying 
on direct tuning of the qubit transition frequency will be 
further discussed elsewhere. 



2-QUBIT GATES: VIRTUAL QUBIT-QUBIT 
INTERACTION 



In this section, we expand the discussions of Ref. [14[ 
on two-qubit gates using virtual excitations of the res- 
onator. To minimize dephasing, we will work with both 
qubits at charge degeneracy (An g = 0). In the ro- 
tating wave approximation, the starting point is there- 
fore Eq. (|2.4|) . To avoid excitation of the resonator, we 
assume that both qubits are strongly detuned from the 
resonator |Aj| = \u aj — u> r \ ^ gj. In this situation, we 
adiabatically eliminate the resonant Jaynes-Cummings 
interaction using the transformation 



U = exp 



Ai 



acr +2 ) 



(5.1) 
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To second order in the small parameters gj/Aj, this 
yields [3, M EH, E3 



H 2q « Wj-o+a + ^ ~2~ az J 



(5.2) a> 



2AiA 2 



(<r +1 cr_ 2 + cr_ 1 £r +3 ) , 



where, as in section mil A[ we have assumed that the cav- 
ity is in the vacuum state and have taken & aj = uj aj + \j ■ 
It is simple to generalize the above expression for an ar- 
bitrary number of qubits coupled to the same mode of 
the resonator. The last term in the above Hamiltonian 
describes swap of the qubit states through virtual inter- 
action with the resonator. Evolution under this Hamilto- 
nian for a time t = 7rAi A2/2gi<72(Ai + A2) will generate 
a V iSWAP gate 14[. This gate, along with the single 



qubit gates discussed in section Hill form a universal set 
for quantum computation [63|. 

In the situation where the qubits are strongly detuned 
from each other, energy conservation suppresses this flip- 
flop interaction. This is most easily seen by going to a 
frame rotating at uj a . for each qubit. In this frame, when 
the qubits are strongly detuned, the interaction term is 
oscillating rapidly and averages out. In this situation, 
the effective qubit-qubit interaction is for all practical 
purposes turned off. On the other hand, for Q ai = uj a2 , 
the interaction term does not average and the interaction 
is effective. 

To turn on and off this virtual interaction, it is neces- 
sary to change the detuning between the qubits. There 
are several ways to do this in the circuit QED architec- 
ture. One possible approach is to directly change the 
transition frequency of the qubits using, as described in 
section [TTl using flux or voltage as control parameters. 
However, as can be seen from eq. (|3.19p . moving the gate 
charge away from the sweet spot will rapidly increase the 
dephasing rate [!, Q . 
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FIG. 5: (Color online) Master equation simulation of the ex- 
cited state population of qubit one (blue) and qubit two (red) 
as a function of time. Full lines: no relaxation and dephasing. 
Dashed lines: finite relaxation and dephasing. The qubits are 
initially detuned form each other with no drive applied. A 
drive, bringing the shifted qubit transition frequencies in res- 
onance, is applied at the time indicated by the vertical dashed 
line. At this time, the two-qubit state start to swap. The pa- 
rameters are: gi/2n = 80 MHz, ff2/27r = 120 MHz, Lu r /2iT = 
5 GHz, Lu ai /2ir = 7.1 GHz, Lu a j2iv = 7.0 GHz, w dl /27r = 
5200.14 MHz. For the dashed lines, cavity decay was taken 
as k/2tt = 0.22 MHz and the qubits decay rates were taken to 
be identical and equal to 71/2-71 = 0.02 MHz and 'y^/2-n = 0.31 
MHz. 



effect of the drive is to shift the qubit transition to 



2A a] 



2 9 -±- 
A' 



(V a) + I) , (5.3) 



where A^- 



,/2A aj — Lu r is the shifted detuning 
also entering in the renormalizcd swap rate. 

The strategy is then to chose u ai 7^ uj a2 such that the 
swap gate is effectively turned off in the idle state. The 
interaction is turned on by choosing a drive amplitude 
and frequency such that ui'^ = u>''. This condition can 
be satisfied with a single drive provided that 31 7^ <?2- A 
master equation simulation of this off-resonant ac-Stark 
tuning is illustrated in Fig. [5] 



The off-resonant ac-Stark shift discussed in sec- 
tion IIII Bl provides another way to tune the qubits in 
and out of resonance. In this situation, one must gener- 
alize the Hamiltonian of Eq. (|5.2[) to include off-resonant 
microwave fields. This is done in appcndix|B]in the pres- 
ence of three independent fields and two qubits. Two of 
the fields, of amplitudes e\ and £2, are used to coher- 
ently control the state of the qubits while the third, of 
amplitude £3 is used to readout the state of the qubits. 

In this section, it is sufficient to take into account a 
single drive £1, assumed to be strongly detuned from 
any resonances. The resulting effective Hamiltonian [5th 
term of Eq. (|B10|) ] contains the swap term already ob- 
tained in the absence of coherent drive in Eq. (|5.2p . The 



B. AC-Dither 

In the same way as the off-resonant ac-Stark shift, ac- 
dither discussed in section HlI Dl can be used to effectively 
tune on resonance a pair of qubits (ac-dither being a low 
frequency version of the off-resonant ac-Stark shift). In 
this situation, the dither frequency must be faster than 
the swap rate between the qubit but still adiabatic with 
respect to the qubit transition frequencies. Moreover, 
similarly to the off- resonant ac-Stark shift, both qubits 
will be blue shifted by the ac-dither. The qubits can 
nevertheless be tuned in resonance by taking advantage 
of their different direct capacitance to the input or output 
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ports (as discussed in section Hi A[) and by using different 
Ec/Ej ratios. 



C. FLICFORQ 

Another approach to tune off-resonant qubits is to use 
the so-calied FLICFORQ protocol (Fixed Linear Cou- 
plings between Fixed Off- Resonant Qubits) [E^. In this 
protocol, one is interested in tuning the effective interac- 
tion between a pair of qubits that arc interacting through 
a fixed linear coupling. The coupling is assumed to be 
off-diagonal in the computational basis such that when 
the qubits are detuned, the coupling is only a small per- 
turbation and can be safely neglected. The interaction 
is effectively turned on by irradiating each qubit at its 
respective transition frequency and choosing the ampli- 
tude of the drives such that one of the Rabi sidebands for 
one qubit is resonant with a Rabi sideband of the other 
qubit. In this situation, the effective coupling becomes 
first order in the bare coupling. 

In circuit QED, FLICFORQ can be used in the disper- 
sive regime to couple a qubit to the resonator or to cou- 
ple a pair of qubits together. As shown in appendix El 
the resonance condition for the first case is A = —Qr 
and this leads to the effective qubit-resonator coupling 
(g/2)(aJ(T- + acr+), at the charge degeneracy point and 
in the RWA. The interaction is first order in the coupling 
g, but of reduced strength. 

Similarly, two qubits that are dispersively coupled to 
the resonator and detuned from each other can be cou- 
pled using FLICFORQ. For simplicity, we again work at 
the charge degeneracy point for both qubits and use the 
rotating wave approximation on the qubit-rcsonator cou- 
plings. To turn on the interaction, two coherent drives 
of frequency Ud- ~ co a . are used. Since the qubits are 
irradiated at their transition frequency, the results of ap- 
pendix [5] cannot be used directly. The corresponding 
effective Hamiltonian is derived in appendix |Dj At one 
of the possible sideband matching conditions and in a 
quadruply rotating frame (see Appendix |D|) . the result- 
ing effective Hamiltonian is 

9l g 2 (A[+A 2 ) 



H 



FF 



uj r a' a+ 



where A' = 



16 a; A' 2 

— oj r with 



(<T yi (T, 



U) 



Ol(2) 



"1(2) 



ft 2 

j -^12(21) 
^ a 12(21) 



■ ), (5-4) 



(5.5) 



the shifted qubit frequency. In this expression, we have 
introduced Qr k — 2<7j£/ c /(uv — ^d k ) the Rabi frequency 
of qubit j with respect to drive k and A ajk = uj aj — Ud k ■ 
This effective qubit-qubit coupling is sufficient, along 
with single qubit gates, for universal quantum computa- 
tion. Moreover, as expected from Ref. [59|, in the FLIC- 
FORQ protocol, the qubit-qubit coupling strength is re- 
duced by a factor of 8 with respect to the bare coupling 
strength. 



D. Fast entanglement at small detunings 

The rates for the two-qubit gates described above are 
proportional to g x (g/A), where g/A must be small for 
the dispersive approximation to be valid. An advantage 
of the dispersive regime is that the resonator is only virtu- 
ally populated and therefore photon loss is not a limiting 
factor. The drawback is that unless g is large, dispersive 
gates can be slow. 

It is interesting to see whether this rate can be in- 
creased by working at smaller detunings. In this situa- 
tion, residual entanglement with the resonator can lead 
to reduced fidelities. As an example, we take for simplic- 



ity 9i = 92 = 9 and w Ql 



UJn 



oj a . Choosing 



2V2g 



1 



(5.6) 



where m > 1 is an integer, one can easily verify that 
starting from \geO) or \egO) at t = 0, the qubits are in a 
maximally entangled state and the resonator in the state 
|0) after a time 



A ' 



(5.7) 



Two-qubit entanglement can therefore be realized in a 
time ~ 1/g and with a small detuning A ~ g, with- 
out suffering from spurious resonator entanglement when 
starting from \ge0) or \eg0). 

It is also simple to verify that \gg0) only picks up a 
phase factor after time T. However, starting from |ee0) 
leads to leakage to \gg2) at time T and therefore to un- 
wanted entanglement with the resonator. As a result, 
while the simple procedure described here is not an uni- 
versal 2-qubit primitive for quantum computation, it can 
nevertheless lead to qubit-qubit entanglement at a rate 
which is close to the cavity coupling rate g. We will de- 
scribe in Section IVII Dl an approach based on Ref. [58[ 
to prevent this type of spurious resonator-qubit entan- 
glement. 

Since the cavity is populated by real rather than vir- 
tual photons during this procedure, it is important to 
estimate photon loss. Starting from \ge0) or \eg0), 
the maximum photon number in the cavity is given by 
"max = (4m 2 — l)/8m 2 . The worst case scenario for 
the rate of photon loss is then T = n max K. The gate 
quality factor, considering only photon loss, is there- 
fore Q ~ 2irT- 1 /T = 32V2gm 2 /[k (4m 2 - l) 3 / 2 ]. For 
g/2n = 200 MHz, k/2tt = 1 MHz and m = 1, this yields 
a large quality factor of Q ~ 1700. 



E. Summary 

The gates discussed in this section (apart from 
Sec. IV Dp rely on virtual population of the resonator. 
A disadvantage of these gates is that they will typically 
be slow. All of them roughly go as g x (g/A) where g/A 
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is a small parameter. Taking g/A ~ 0.1 and g/2ir = 200 
MHz, we see that the rates of these gates will realisti- 
cally not exceed 20 MHz. Although not very large, this 
rate nevertheless exceeds substantially t he ty pical decay 
rates 71 , 7^ and k of circuit QED 0, [15|, ll6| and should 
be sufficiently large for the experimental realization of 
these ideas. An advantage of these virtual gates is how- 
ever that, since the resonator field is only virtually pop- 
ulated, the gates do not suffer from photon loss. As a 
result, these gates could still be realized with a resonator 
of moderate Q factor (which is advantageous for fast mea- 
surement 0). 

It is also interesting to point out that, in the situation 
where there are more than two qubits fabricated in the 
resonator, the same approach can be used to entangle 
simultaneously two or more pairs of qubits. This is done, 
for example, by taking uj ai = u a2 7^ ui a3 = uj ai while still 
in the dispersive regime. It is simple to verify that this 
corresponds to two entangling gates acting in parallel on 
the two pairs of qubits. This type of classical parallelism 
is an important requirement for a fault-tolerant threshold 
to exist [H[. 

Finally, we point out that the dispersive coupling can 
also be used to couple n > 2 qubits simultaneously. This 
is done by tuning the n qubits in resonance with each 
other but all still in the dispersive regime. This leads to 
multi-qubit entanglement in a single step. 



VI. 



CONDITIONAL ENTANGLEMENT BY 
MEASUREMENT 



As discussed in section MI CI and in more detail in 
Ref. [1, [3, [HI, [H, [l|| , measurement can be realized in 
this system by taking advantage of the qubit-state de- 
pendent resonator frequency pull. In the presence of 
a single qubit, the resonator pull is x a z and becomes 
Y^Jj Xj a zj in the n-qubit case. If the pulls Xj are dif- 
ferent for all j's and large enough with respect to the 
resonator linewidth k, each of the different n-qubit states 

{\999-9), \egg--g), \g e 99--g), ■■■> \eee...e)} pulls the cav- 
ity frequency by a different amount. In this situation, it 
should be possible to realize single shot QND measure- 
ments of the n-qubit state. In a test of Bell inequalities, 
this multi-qubit readout capability would offer a power- 
ful advantage over separate measurement of each qubit. 
Indeed, in the latter case the effective readout fidelity 
would be the product of the individual readout fidelities. 

The situation is also interesting in the case where the 
pulls x are equal. For example, in the two qubit case, 
when xi = Xi = X (and taking u ai = uj a2 = to a for 
simplicity), the dispersive Hamiltonian of equation (|5.2p 
can be written as 

H 2q « [ui r + xfazi + <Tz 3 )] a^a + ^ - (w a + x) a Xj 



j=l,2 



While this Hamiltonian is not QND for measurement of 
a Zj , it is QND for measurement of (a Zl +c Z2 ) since [(cr Zl + 
<Jz 2 ),H 2q \ = 0. 

More interestingly, in this situation, the states 
{|<?e), |eg)} while they may have different Lamb shifts 
have the same cavity pull. This implies that they cannot 
be distinguished by this measurement. The consequence 
of this observation can be made more explicit by rewrit- 
ing (|6.ip as 



H 2q = UrcJa + x i\ip+) (Vm 



2 X ( a f a + - 



(|ee)(ee| - \gg)(gg\), 
(6.2) 



where \ip±) = (\ge) ± \eg)) /y/2 are Bell states. As a re- 
sult, the projection operators corresponding to measure- 
ment of the cavity pull are Tl g = \gg)(gg\, n e = |ee)(ee| 
andn Bc n = + with Y, k=g e Bdi u k = 

1. 

An initial state will thus, with probability 

(-0i|nBcii|'0i): collapse to a state of the form \t/jf) = 
c+\ip+) + C-\ip-}. Since the Bell states \ip±) are eigen- 
states of H 2q , further evolution will keep the projected 
state in this subspace. For certain unentangled initial 
states [e.g. (\g) + |e)) ® (\g) + |e))/2, which is created by 
7r/2 pulses on each qubit], the state after measurement 
will be a maximally entangled state. As a result, condi- 
tioned on the measurement of zero cavity pull, Bell states 
can be prepared. This type of entanglement by measure- 
ment for solid-state qubits was also discussed by Ruskov 
and Korotkov (65( . We also point out that entanglement 



+ X(cr +1 CT_ 2 +cr_ 1 cr +2 ) 



(6.1) 



by measurement and feedback was studied in Ref [66 



VII. 2-QUBIT GATES: QUBIT-QUBIT 
INTERACTION MEDIATED BY RESONATOR 
EXCITATIONS 

In this section, we consider gates that actively use the 
resonator as a means to transfer information between the 
qubits and to entangle them. More precisely, we will take 
advantage of the so-called red and blue sideband transi- 
tions. We first start by a very brief overview of ion-trap 
quantum computing to show the similarities and differ- 
ences with circuit QED and then present various proto- 
cols adapted to circuit QED. We discuss the realization 
of quantum gates based on these sideband transitions in 
section lyTTDl 



A. Ion-trap quantum computation 

Sideband transitions are used very successfully in ion- 
trap quantum computation [6^, [68[. This approach to 
ion-trap quantum computing was first discussed by Cirac 
and Zoller (57j and, more recently, was adapted to two- 
level atoms by Childs and Chuang . Following closely 
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Ref. [HI the Hamiltonian describing a trapped ion is 

H = oj r a)a + ~^ a z — \ cos{ri(aJ + a) — ujt]a x , (7.1) 

where uj r is the frequency of the relevant mode of oscil- 
lation of the ion in the trap, u) is the laser frequency, 
r) = \/k 2 /2Nniuj r is the Lamb-Dickc parameter and A 
is the amplitude of the magnetic field produced by the 
driving laser (58j . Assuming 77 1 and choosing the 
laser frequency u = u a + nu) r , the only part of the inter- 
action term that is not rapidly oscillating and thus does 
not average out is 

\n n 

ffn-blue « ( at " CT + + a " ff - ) ■ ( 7 ' 2 ) 

For uj = u> a — ntu r , we get 

i?„-red « (a f 'V_ + a n a+) . (7.3) 

These correspond to blue and red n-phonon sideband 
transitions, respectively It is known that these tran- 
sitions, along with single qubit gates, are universal for 
quantum computation on a chain of ions [58| . 

We first note that the rate of the sideband transitions 
scales with A which itself depends on the (variable) laser 
amplitude. In circuit QED, we will see that the (fixed) 
qubit-rcsonator coupling g takes the role of the param- 
eter A. Moreover, while the Hamiltonian (|7.ip allows 
for all orders of sideband transitions, we will see that 
when working at the charge degeneracy point, the sym- 
metry of the circuit QED Hamiltonian only allows side- 
band transitions of even orders. Finally, it is interesting 
to point out that the Hamiltonian (|7.ip also describes 
a superconducting qubit magnetically coupled to a res- 
onator. In that situation, the qubit would be fabricated 
at an anti-node of current in the resonator. The zero- 
point fluctuations of the current generate a field that 
couples to the qubit loop. For example, a superconduct- 
ing charge qubit magnetically coupled to the transmis- 
sion line resonator would have in its Hamiltonian a term 
of the form A cos [r](aJ + a) + 0\ a x , where A = Ej and 
7] = 7rM/° ms /$o with M the mutual inductance between 
the center line of the transmission line, 7° ms the rms value 
of the vacuum fluctuation current on the center line in 
the ground state [H, [6^] , $0 is the superconducting flux 
quantum and 9 is controlled by the dc flux bias. 

An important difference between trapped ions and the 
superconducting flux qubit analog is that in the first case 
A scales with the laser amplitude while, in the second 
case, it is equal to the Josephson energy. The latter can 
only be so large in practice and is difficult to tune rapidly. 
Finally, we point out that solid-state qubits, in contrast 
to ions, have a natural spread in transition frequencies. 
This allows to address the qubits individually with global 
pulses, without requiring individual bias lines to tune 
them. 



B. Sideband transitions in Circuit QED 

Our starting point to study sideband transitions in 
circuit QED is the two-qubit Hamiltonian of Eq. (|2.3[) . 
Here, we keep the gate charge dependence and do not 
initially make the rotating wave approximation. The ef- 
fective Hamiltonian describing the sideband transitions 
in circuit QED is obtained in appendix [Bl As discussed 
above, in this appendix we consider the presence of two 
qubits and three coherent drives. Two of these drives, of 
frequency u>d 1 2 and amplitude £1,2, are used to drive the 
sideband transitions. The third drive, of frequency u>d 3 
and amplitude £3, is used to measure the state of the sys- 
tem. For simplicity of notation, in this section we focus 
on a single qubit coupled to the resonator and drop the 
j index. 

The red and blue sideband transitions, illustrated in 
Fig. [5J are given by the last two lines of Eq. (|B10[) . These 
correspond to single and two-photon sideband transi- 
tions, respectively Higher order transitions are neglected 
due to their small amplitudes. We rewrite these terms 
here in a more explicit form. This is done by going to a 
frame rotating at uj r for the resonator and at the shifted 
frequency w'^ for the qubit, where 

<=-+3t+S: + 4(<^)- (7 - 4) 

Here we keep the contribution of {a* a) to the frequency 
shift in order to take into account the presence of the 
measurement drive £3. When evaluating (a^a), it is im- 
portant to remember that, a'W is in the displaced frame 
defined in appendix [Bl Setting Wd ,,, = oj'' + ui r the only 

k { ' 

term that does not average to zero due to rapid oscilla- 
tions is the third line of Eq. (|B10[) which gives 

Hn = -C3^ (0 + a + a_ at) , (7.5) 

where k = 1 or 2. Following the notation introduced 
in Sec. Ill A[ we use c = cos 9 and s = sin# with 9 the 
mixing angle. The above Hamiltonian corresponds to 
the one-photon red sideband transition. Alternatively, 
taking w<j , = cj" — u r , we obtain the one-photon blue 
sideband transition 

flbi = - c 5^ ( CT + at +*■-*)■ (7.6) 

On the other hand, the last line of (|B10|) will not aver- 
age to zero if the drive frequencies are chosen such that 

uj dk ± uJ dk , = lu" + UJ r (7.7) 

or 

2u dhil) =u%+u> r . (7.8) 
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FIG. 6: (Color online) Red and blue sideband transitions, a) 
One-photon transitions, b) Two-photon transitions. 

With these choices of drive frequencies, we obtain 

O o (7 ' 9) 

#b2 = s 3 o~A OA ' CT+a + f7 - a ) • 

corresponding to two-photon red and blue sideband tran- 
sitions, respectively. These two-photon transitions are 
illustrated in Fig. [5b). 

Because of the dependence of Eqs. (|7.5[) and (|7.6[) on 
the cosine of the mixing angle 9, the first order red and 
blue sideband transitions are forbidden at the charge de- 
generacy point. As discussed in appendix [El this can be 
linked to the symmetry of the Jaynes-Cummings Hamil- 
tonian. Since it is more advantageous to work at the 
sweet spot to minimize dephasing, these first order tran- 
sitions therefore appear to be of limited interest for co- 
herent control in circuit QED. A similar selection rule, 
for flux qubits coupled to a LC oscillator, was noted in 
Ref. [2(j. There, it was suggested to work with single- 
photon sidebands by biasing the flux qubit close to de- 
generacy, but not exactly at the sweet spot. Moreover, 
selection rules for a flux qubit irradiated with classical 
microwave signal were also studied in Ref. [7(| . We also 
point out that, since the frequencies corresponding to 
these first order transitions (w", ± uj r ) are in practice 
very detuned from the resonator frequency tu r , signals at 
these frequencies would be mostly reflected at the input 
port of the resonator. Very large input powers would 
therefore be require to compensate the attenuation. 

On the other hand, because of their sin 8 dependance, 
the two-photon transitions have maximal amplitude at 
the sweet spot. Moreover, since in this case we require 
the sum or the difference of the drive frequencies to match 
the sideband conditions, these frequencies can be chosen 
such that there is minimal attenuation (while still avoid- 
ing measurement- induced dephasing [3 Gil Ell ) ■ How- 
ever, the rate for these 2-photon transitions is of order 
g x (£ln/2A a ) 2 , which is g times the square of a small pa- 
rameter. As is discussed in section IV El obtaining large 
rates will require large coupling strength g. 

To realize logical gates based on the red and blue side- 
bands, these transitions must be coherently driven. The 
corresponding simulated coherent oscillations arc illus- 
trated in Figurc[7]for the two-photon red (a) and blue (b) 
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FIG. 7: (Color online) (a) Red and (b) blue coherent side- 
band oscillations. The system parameters are uj a /2n = 7 
GHz, u r /2-K = 5 GHz, g/2n= 100 MHz. For the color curves, 
the damping is k/2tt = 0.22 MHz, 7i/2tt = 0.02 MHz and 
■7^ /2-7T = 0.31 MHz. These last two values are taken from 
the measured value reported in Ref. Q. For the gray curves, 
damping was set to zero. The value of the drive frequencies 
and amplitudes are: a) ujdi/2n = 4593.59 MHz, ujd 2 /2n = 
6650.605 MHz, £i/2tt = 841.91 MHz and £ 2 /2tt = 843.417 
MHz, b) uj dl /2ir = 4598.39 MHz, uj d j2n = 7476.47 MHz, 
£i/2tt = 1241.51 MHz and £ 2 /2tt = 1252.56 MHz. 

sideband transitions. In these numerical calculations, the 
drive frequency ui^ is chosen ~ 400 MHz away from the 
resonator frequency to avoid measurement-induced de- 
phasing of the qubit. In a first step, the second drive fre- 
quency ujd 2 is then chosen u sing the condition of Eq. (|7.7p . 
Using simulated annealing [7l| , the drive frequencies and 
the corresponding drive amplitudes £1(2) ar e then varied 
to optimize the fidelity of the sideband transitions. The 
best (but not necessarily optimal) values obtained are 
given in the caption of Fig. [71 Using these parameters, 
we obtain a population transfer of 0.83 for the red side- 
band and of 0.86 for the blue sideband (without damping, 
we obtain near perfect population transfer of 1.0 in both 
cases). This relatively low population transfer is essen- 
tially due to the small dephasing time 1 /j2 = 500 ns used 
here with respect to the slow rate ~ 5 MHz of the red and 
blue sideband transitions. It is interesting to point out 
that this value of the rate is about 4 times bigger than 
expected from the perturbative estimates of Eq. (|7.9p . 
These estimates should be taken as lower bounds which 
can be improved by numerical optimisation. 

C. One-photon sidebands at the sweet spot using 
ac-dither 

In the previous section, we saw that the symmetry of 
the circuit QED Hamiltonian does not allow for first or- 
der red and blue sideband transitions at the charge de- 
generacy point. However, with ac-dither discussed in sec- 
tion IIII D[ it is possible to take advantage of the small 
gate charge excursions away from the sweet spot to ob- 
tain a one-photon transition while staying, on average, 
at the sweet spot. This can be seen as a low frequency 
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version of the two-photon transitions. 

To analyse this situation, we focus on a single qubit 
coupled to the resonator and in the presence of a single 
coherent drive of amplitude e and frequency Ud ■ At the 
charge degeneracy point (rig C = 1/2) and with ac-dithcr 
on the voltage port, we get for the Hamiltonian: 



H = oj r a) c 



E 3 A , , 
a z + — cos(ui ac t)a x 



2 2 

+ g (2rig c cos(w ac f) + a x ) (a) + a) 

+ e{a*e- iWdt + ae +luJdt ), 



(7.10) 



where we have defined A = 8Ecrig C and have taken the 
ac-gate bias to be rig C cos(cj ac f). Assuming that the dither 
frequency is small with respect to the qubit transition 
frequency but large compared to the coupling strength, 
g <C c<j ac <C uj a , we move to the adiabatic basis to obtain 

H = w r <Ja + - jl y 1 &z + e(ah- iUdt + ae +i " dt ) 

+ g [2rig C cos(w ac f) + sin 0(f)<r z + cos Q(t)a x ] (a* + a) , 

(7.11) 

1 /2 

where uj^ d (t) = [Ej + A 2 cos 2 (w ac f)] is the instanta- 
neous qubit transition frequency. 

Comparing with Eq. (|2.3j) , we have the following map- 
ping 

[i i— > — 2rig C cos(w ac f), ci- > sin 0(f), si— >— cos 0(f) 

(7.12) 

and it is possible to use directly all of the results obtained 
in the previous section. The important point is that we 
now get a b z component even at the charge degeneracy 
point, which means that the red and blue sidebands will 
be allowed to first order. This is due do the small excur- 
sions away from degeneracy provided by the ac-dithcr. 

Assuming that the dither amplitude is small with re- 
spect to the bare qubit transition frequency Ej , we take 
0(f) ?s Acos(cu & ct)/ Ej which will yield simple Bcsscl 
function FM modulation sidebands for the qubit tran- 
sition. Following section IIII Dl and using the results of 
the previous section, we get the blue sideband transition 
for uj d = u'l + uj r ± w ac 

/ n R \ T {8E c n a a 



Hbi B 



9 



\2A a 



Ji 



E, 



a) ij + + ao"_ 



(7.13) 



and for lu ( i 
sition 



uj r ± aj ac the red blue sideband tran- 



a f CT_) . (7.14) 



Here we have focused on the first dither sideband J\ (z) 
only. With respect to the two-photon transitions of 
Eq. (|7.9[) . we have simply replaced a factor of f2fl/2A 
by the ac-sideband modulation J 1 (8Ecrig C /Ej). Both of 
these quantities are smaller than unity, so whether this is 
advantageous depends on the parameters of the system. 
As an example, taking Eq ~ 5 GHz, Ej ~ 6 GHz [l6| . 
and a 10% of 2e excursion for the dither, n a g c — 0.1, we 
have J x {8E c rC/E 3 ) - 3/10. 



D. CNOT from sideband transitions 

In this section, we show how to obtain non-trivial two- 
qubit gates from red and blue sideband transitions. This 
is included for completeness, with most of the results 
already known from Cirac and Zollcr [571 for three level 
atoms and from Childs and Chuang [58[ for two level 
atoms. Here, we follow closely the results and notation 
ofRef. [13 . 

Following Ref. [58|, we introduce the unitary opera- 
tors [zl 



exp 
exp 



z 



, (7.15) 
. (7.16) 



R 



corresponds to a pulse on the blue-sideband for 



qubit j and R 



the red-sideband. Here, d> is the 



phase of the driving field. 

In addition to the above resonator-qubit operations, we 
introduce the singlc-qubit unitary operators correspond- 
ing to the effective Hamiltonians discussed in section Hill 
We denote the single qubit flip (x) and phase (z) opera- 
tors acting on qubit j as: 



R x A0) = exp 



-l— Or, 

2 



Rf(0)=exp 



2 Zl 



(7.17) 

Another useful single qubit unitary is the Hadamard 
transformation (in the basis {|e), \g)}) 



Wl 1 
V2V 1 



(7.18) 



which can be obtained by a rotation at an angle between 
the x and z axis or, equivalcntly, by a sequence of one- 
qubit gates 



= exp 



7T 

—i — 
2 



V2 

Rf( 7 r/2)R|( 7 r/2)Rf( 7 r/2) 



(7.19) 



Before building a universal two-qubit gate from these 
elementary operations, we first discuss a simpler proto- 
col to create conditional entanglement. This protocol is 
based on Ref. [33| and relies on entangling one qubit to 
the resonator and then transferring the entanglement to 
qubit-qubit correlations. This is realized by the sequence 



R-(7r,0)R-(V2,0)Rr(7r,. 



(7.20) 



where the phase <p is arbitrary. It is simple to verify that 
this sequence will leave \gg0) unchanged but will create 
maximally entangled states when starting from \ge0) or 
\eg0). However, the state |ee0) will irreversibly leak out 
of the {|0), |1)} photon subspace and the qubits will get 
entangled with the resonator at the end of the pulse se- 
quence. Hence, while (|7.20[) is not an universal two-qubit 
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primitive, it can nevertheless be used to generate condi- 
tional entanglement. This sequence can also be realized 
with all blue sideband pulses. In this case, |eeO) is left 
unchanged while starting with \ggO) will create spurious 
entanglement with the resonator. 

In the above sequence, the spurious entanglement oc- 
curs in the second step where the initial state |eeO) picks 
up a contribution from the two-photon state \gg2). For 
this state, the evolution is "faster" by a factor of \[2 from 
evolution in the one-photon subspacc. Because of this, 
the last step in the sequence cannot completely undo the 
qubit-resonator entanglement and leaves them partially 
entangled. To solve this problem, Childs and Chuang [58[ 
introduced the qubit-resonator gate 

Pj = R+(-7r/2,0)Rt(7rV^,-7r/2)R+(-7r/2,0). (7.21) 

In the basis { | eO) , |el), \g0), |<?1)}, Pj takes the form 

Pj = diag(l,e-"/ v ^,e l7r/%/5 , -1). This gate entangles 
the qubit with the resonator (because of the minus sign 
in the last element) but it does not lead to leakage 
into higher photon states. Using this gate, Childs and 
Chuang [58| proposed a sequence of red, blue and single- 
qubit operations that generates a CNOT gate. 

Here, building on Eq. (|7.20p . we note a simpler entan- 
gling two-qubit gate (in the basis {|ee), \eg), \ge), | <?<?}}) 



IV = Rf(7r,^P 2 R+(7r, 



= diag(l,e 



:tt/V2 -i 

1 t- 



7V2 



(7.22) 



Using this gate, it is possible to obtain a CNOT gate 
which relies only on single qubit unitaries and blue- 
sideband transitions: 

CNOT = Rf([l + v^]7r)H 2 U v ,l^(-7r/v^)H2. (7.23) 

Because it relies on Pj, this gate does not lead to un- 
wanted qubit-resonator entanglement. 



E. Summary 

The gates presented in this section rely on real excita- 
tions of the resonator to mediate entanglement between 
the qubits. These gates are based on perturbation theory 
and are therefore relatively slow. For example, for the red 
and blue sideband oscillations studied in section IVII Bl 
we have found after numerical optimization rates of ~ 
5 MHz. With a dephasing time of 500 ns, this corre- 
sponds to a quality factor of about 9, larger than what 
was expected from perturbation theory. While this qual- 
ity factor is not large enough for large scale quantum 
computation, it is certainly enough to demonstrate the 
concept experimentally. 

Finally, a disadvantage of the gates based on real ex- 
citation of the resonator is that they are susceptible to 
photon loss and therefore require relatively large Q res- 
onators. This conflicts with the requirement of fast read- 
out. 



VIII. GEOMETRIC PHASE GATE 

The gates discussed in the previous sections were based 
on real or virtual transitions. In the present section, 
we discuss a different approach, based on the geometric 
phase. This was already discussed in the context of ion- 
trap quantum computing [46l \4l\ . This gate is based on 
the first term of the second line of the Hamiltonian (|B10[) : 



H = - 9jBjCT Zj (a f + a), 



(8.1) 



where Bj is given by Eq. (|B9[) . Here, we work at the 
charge degeneracy point where C1.2 = 0. Although the 
Hamiltonian (|8.1[) docs not couple the qubits directly, it 
couples both qubits to the resonator field a. By using 
a time dependent drive on the resonator, and hence dis- 
placing the field a in a controlled manner, it is possible 
to induce indirect qubit-qubit coupling without residual 
entanglement with the field. 

To see this explicitly, we first rederive the effective 
Hamiltonian (|8.1[) in the presence of a single drive of fre- 
quency to and of amplitude e(t). Here, we allow the am- 
plitude to be time-dependent and complex. Moreover, we 
will assume that the qubits are detuned from each other 
(w^ ^ w„ ) such that the flip-flop interaction can be 
neglected and that both qubits are dispersively coupled 
to the resonator. In a frame rotating at the resonator 
frequency u> r , we obtain 



#(*) « \ E (< (*) + 2 /^) at + 2 /; (*)< 



J'=1.2 



(8.2) 



where 



gh~ lA ^ 



Jo 



(8.3) 



Because of the time dependent drive, the shifted qubit 
transition frequencies are now time-dependent. 

Note that, following the procedure of appendix [5] we 
have added the effect of cavity damping k to Eq. (|8.3I) . 
Since 



[H(h),H(t 2 )] 

= 2i (lm[f*(t 1 )f 2 (t 2 )}+lm[f;(t 1 )f 1 (t 2 )})a z 
= 2iF(t 1 ,t 2 )a Zl a Z2 



a Zo 



(8.4) 



commutes with H(t) for all times, evolution under the 
effective Hamiltonian (|8.2|) is given by 

U(T,0) = e -iSo ^^e-kSo Sodtdt'[H(t),H{t')]_ ( 8 .5) 

To avoid unwanted entanglement of the qubits with the 
resonator at the final gate time T, we choose the time- 
dependent drive amplitude e(t) such that [-ill |47| 



dtf j (t)=0. 



(8.6) 
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With this choice of e(t), evolution under H(t) corre- 
sponds to the application of a local phase shift 

4>j{T) = f T dtuj'^t) (8.7) 
Jo 

to each qubit and of a conditional phase shift 
exp [-i$ 12 (t)a zl cr Z2 /2} where 

$ 12 (T) = 2 [ [ dtdt'F(t,t'). (8.8) 
Jo Jo 

For < &i2(r) = ±7r/2, this two-qubit operation is known to 
be equivalent to the CNOT gate, up to one qubit gates. 

Our goal is therefore to choose the drive e{t) such that 
the qubits accumulate a phase $ 12 (T) = ±7r/2 in the 
smallest time T possible. This minimization has to take 
into account several constraints. In addition to (|8.6p , we 
take the drive to be off at the start and at the end of the 
gate: 

e(0)=e(T) = 0. (8.9) 

The assumption that the drive is initially turned off is 
already built into our expression for fj (t) in Eq. (|8.3[) . To 
prevent further phase accumulation after the gate time 
T, we also require that 

a(T) = 0, (8.10) 
d(T) = 0, (8.11) 

where a(t) is the amplitude of the classical field inside 
the resonator and is given by Eq. (|3.4|) [or Eq. (|A5|) in 
presence of cavity damping]. From Eq. (|3.4p . we have 
that a(T) = will be automatically satisfied if Eqs. (|8.9p 
and (|8.10p are satisfied. 

While satisfying the above constraints, the external 
field e(t) must also be chosen such that the approxima- 
tion that led to Eqs. (|8.2p and (|8.3p are valid. There are 



two approximations. The first one, as in Section UlI B[ is 
that we adiabatically eliminated transitions in the qubit 
caused by the external drive. This requires the drive 
to be either of small amplitude or sufficiently detuned 
from the qubits. The second approximation is the dis- 
persive transformation, which here only has the effect of 
renormalizing the qubit transition frequency lu".. This 
approximation breaks down on a scale given by the crit- 
ical photon number n cr it = A 2 /4g 2 [lj, 0- We must 
therefore choose the drive amplitude and frequency such 
that the resonator never gets populated by more than 
?T-crit photons. Beyond this number, additional mixing 
between the qubits and resonator states is possible and 
would likely lead to unwanted qubit-resonator entangle- 
ment [7^] 

Finally, we also require that the measurement-induced 
dephasing time of the qubits due to the control drive, 
the gate-induced dephasing time, to be smaller than 
the 'intrinsic' dephasing time T 2 . An expression for 
the measurement-induced dephasing rate is give n in 
Eq. (j3~T2)) . Here, we will use Eq. (5.16) of Rcf. [3 that 
is more appropriate for a time-dependent drive. As was 
discussed in section IIII CI and in Ref. [l8[ , by working 
at sufficiently large detuning A r , gate- induced dephas- 
ing can be made negligible while still maintaining a large 
rate for the gates. 

It is useful to take into account these constraints by 
developing the pulse envelope e(t) in Fourier components 



oo 

e(*) = c ^ +luJnt - (8-12) 



Using this expression, we rewrite fj(t) as 



n— — oo 
oo 

~ ~Xj E/ c " 



-Kt/2 



-iAjt 



(Aj + A r + cj„)(A r + uj n - in/2) (Aj + i«/2)(A r + u> n - in/2) (A 3 + A r + (J n ){Aj + in/2) 

e i(A r +o!„)t _ g-ret/2 



(A.,. + u) n — in/2) 



where we have assumed the resonator-qubit detuning Aj 



is large, such that Aj 3> {A,. 



k/2}, and have 



dropped the small and fast oscillating last term in the 
second expression. This approximate expression will be 
useful in obtaining analytical estimates for the expected 
gate time T. 



.13) 



r 



Using the approximate expression for fj (t) , we can now 
recast the above constraints in more simple forms. Us- 
ing (|8.13p . the no-effect on the resonator constraint of 
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Eq. (|8.6[) can be written as 



(i 



dtf 3 (t) 



-xj Cn 

n— — oo 
oo 

Xj ^ ] c n A r 



e !(A r +u„)t 



o (A, 



UJ r , 



in/2) 



.14) 



Without the approximation for /j(i), the coefficient A„ 
depends on the qubit index j which means that there 
would be one extra constraint to satisfy. 

Moreover, to satisfy Eq. (|8.9p . we take uj n = 2im/T 
which implies that 



e(0)=e(T)= J2 c « = - 



.15) 



To satisfy Eq. ([8TT0]l . we have 

e i(A r +^„)T _ e -KT/2 



(A r + oj n — in/2) 

oo 

CnB n =0. 



(8.16) 



Casting, in a simple form, the constraints that there 
is not qubit transition and that the dispersive approxi- 
mation holds is more difficult. For the former case, we 
require that 



\f j (t)e- iA r t /g j 



AjA r 



(8.17) 



be smaller than about 1/10 for all time. In terms of the 
drive amplitude we thus require 



\e(t)\ 



< 



10 9j 



(8.18) 



Finally, for the dispersive approximation to hold, we re- 
quire that the average photon number populating the 
resonator be no larger than the critical photon number: 



A 1 



or, equivalently, 



\e(t)\ 



< 



AjA r 



2g 3 



.19) 



.20) 



In practice, we therefore only have to deal with Eq. (|8.18p 
and only with the qubit for which \ Aj/gj \ is smallest. In 
the dispersive regime and in the k = approximation, 
we thus require \e(t)\ < A r which means that there is no 
gate effect at A r — > 0. 

Of all the constraints, the last one is the most diffi- 
cult to deal with because it must hold at all intermediate 




~ a 




10 20 30 

A r /2ji [MHz] 



20 40 60 80 
A r /2jt [MHz] 



FIG. 8: (Color online) a) Gate time (full blue and green 
dashed) and gate-induced dephasing (red dotted) as a func- 
tion of detuning A r . The green dashed line is obtained ana- 
lytically by decomposing the pulse e(t) over four Fourier co- 
efficients and using the constraints derived from the approx- 
imate expression in Eq. (|8.13|l . The full blue line is obtained 
numerically by using the constraints derived from the exact 
expression in Eq. (|8,13[) and decomposing over 5 Fourier co- 
efficients. The dotted red line is the gate-induced dephasing 
time as calculated from Ref. and presented on a log scale. 
The system parameters used are given in the text, b) Quality 
factor of the geometric phase gate as a function of detuning. 
The red dashed line takes into account gate-induced dephas- 
ing and an additional T2 of 500 ns. For the full blue line, we 
have taken I/T2 = 0. Clearly, gate-induced dephasing is not 
a limitation in practice. 



times t and because it involves the absolute value of the 
complex drive amplitude e(t). It is nevertheless possible 
to find an analytical expression for the gate time T by 
developing on four Fourier components {c_x, Co, ci, C2}. 
In this situation, all constraints can be solved for ana- 
lytically and it is possible to find an expression for the 
gate time. This expression is however unsightly and we 
only give here relevant limits. In the situation where A r 
is large, we find that 



T 



TTjArl 
XlX2 



.21) 



In the dispersive approximation this is roughly given by 
T - 1007r|A r |/.9 2 . In the small A r limit we rather find 
that T (x 1/1 A,. I . As explained above, this behavior is 
expected from Eq. (|8.18[) which is saying that the field 
amplitude goes to zero at zero detuning in the k = 
limit. The full expression for T as a function of detuning 
A r is plotted in Fig. [5^,) [dashed green line] . The system 
parameters used here are are gi/2n = .92/2^ = 100 MHz, 
Ai/2tt = 1 GHz and A 2 /2tt = 1.1 GHz. The approxi- 
mate expression for T at large A r is of similar form to 
that already obtained for gates based on the dispersive 
approximation: a rate equal approximately to the square 
of the coupling g over a detuning. Here however the de- 
tuning is A r and not the qubit-resonator detuning A. In 
the dispersive approximation, the former can however be 
made much smaller than the latter and this gate could 
in principle be faster than the dispersive based gates dis- 
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cussed in the previous sections. 

Going beyond the assumptions made to obtain the ap- 
proximate expression for fj(t) in Eq. (|8.13p . we solved 
numerically the optimization problem. Without this ap- 
proximation there are now five constraints to optimize 
over [the constraint Eq. (|8.14|) now has to be satisfied in- 
dependently for both qubits] and thus a minimum of five 
Fourier coefficients have to be used. The full blue line in 
Fig. [gji) shows the numerically found gate time as a func- 
tion of detuning and using the same system parameters 
as given above. As can be seen from this figure, going 
beyond the approximations used to get an analytical es- 
timate and increasing the number of Fourier components 
improved significantly the gate time. In this situation, we 
get an optimal gate time of T ~ 50 ns. Further optimi- 
sation can be made by increasing the coupling strength, 
and could likely be made by increasing the number of 
Fourier components. 

Figure [8^1) also shows the gate- induced dephasing time 
[dotted red line] associated with the pulse required to 
implement the geometric phase gate. This dephasing 
time is obtained from Eq. (5.16) of Ref. [l|[ and as- 
suming a resonator damping rate of k/2-k — 0.1 MHz. 
At detunings larger than about A r /27r ~ 25 MHz, 
the measurement-induced dephasing time is significantly 
larger than the 'intrinsic' T2 of 500 ns already measured 
111 circuit QED @ and can thus be ignored. However, 
this induced dephasing time goes down rapidly with de- 
tuning such that there is a detuning (A r /27r ~ 10 MHz 
with the chosen parameters) below which it is smaller 
than the gate time, meaning that the geometric phase 
gate cannot be used. 

Figure [5Jd) shows the quality factor, as defined by 
Eq. (|3.14j) and where the dephasing rate was taken to 
be the sum of the gate-induced dephasing rate and of the 
'intrinsic' dephasing rate I/T2. For the red dashed line, 
we have taken T2 = 500 ns, while T2 was set to infinity for 
the full blue line. For T2 = 500 ns and the present system 
parameters, the quality factor is maximum for detunings 
of about 30 MHz. However, as is clear from the full blue 
line, gate-induced dephasing time is not a limitation in 
practice and the quality factor could be significantly bet- 
ter once T2 is improved in this system and by working 
at larger detunings. Moreover, the actual magnitude of 
the quality factor could likely be improved by increasing 
the number of Fourier components. We also point out 
that the present results have been obtained within the 
approximations used to derive the effective Hamiltonian 
of Eq. (|8.2[) . Significant improvements could be made by 
numerically optimizing the full system's master equation. 
In this case, the results would not be limited by the dis- 
persive approximation used here and we expect a much 
better quality factor. 

Finally, it is important to realize that whenever two 
qubits are present in the resonator and the system is be- 
ing driven, this geometric gate is in action. This could 
lead to unwanted qubit-qubit entanglement and, since 
Eq. (|8.6p might not be satisfied, qubit-resonator entan- 



glement if the drive amplitude and frequency are not cho- 
sen appropriately. 



IX. CONCLUSION 

We have explored the realization of single and two- 
qubit gates in circuit QED, using realistic systems pa- 
rameters. We have shown how all single-qubit rotations 
can be realised with minimal measurement-induced de- 
phasing of the qubit. In this context, two approaches 
to change the qubit transition frequency without dc bias 
away from the sweet spot were discussed. Both rely on 
off- resonant irradiation of the qubit. Interestingly, the 
ac-dither approach could help in reducing dephasing by 
protecting the qubit from low-frequency noise. 

Five types of two-qubit gates were discussed. The 
fastest gate discussed in this paper is based on direct 
coupling of qubits with the resonator. In this paper how- 
ever, we have focused on gates requiring no dc-bias away 
from the sweet spot or additional control lines not present 
in the current circuit QED architecture. These gates can 
in practice be slower but could be implemented in the 
current circuit QED design, without additional design 
elements. Moreover, these gates have the important ad- 
vantage that they do not cause extra dephasing of the 
qubits by moving away from the sweet spot. 

The rates and quality factors for these gates, as ob- 
tained from perturbation theory, are summarized in Ta- 
ble [B These are given for g/2n = 100 MHz and n/2n 
= 0.1 MHz ('current') and g/2n = 200 MHz and n/2n 
= 0.01 MHz ('better'). All other parameters are equal 
and discussed in the table caption. These parameters 
correspond to already realised values and to realistic val- 
ues that could be obtained in future realizations of cir- 
cuit QED, respectively. The quality factors are calcu- 
lated in the same way as in Eq. (|3.14p : the rates divided 
by 2T 2 _1 (the full width at half-maximum of the qubit 
spectral line). However, for the direct coupling and the 
sidebands, the loss of a single photon has a large impact 
and the quality factors are therefore given by the rates 
divided by the mean of the two contributing decay chan- 
nels: (k + 2T 2 ~ 1 )/2. For the geometric phase gate, the 
'current' quality factor, and the corresponding rate, are 
taken from the inset of Fig. [5J The 'better' results are 
obtained from a similar calculation. 

Moreover, we note that the same values of T2 and k 
were used for the direct coupling gate requiring dc excur- 
sions away from the sweet spot and the other gates that 
do not require such excusions. As discussed before, it is 
likely that T2 or k would be reduced due to the dc-bias.. 
The rates and quality factors given on the first line of the 
table can therefore be taken as upper bounds. 

It is also important to realise that, appart from the first 
line of the table, these results have been obtained by per- 
turbation theory and are thus lower bounds on the rates 
and quality factors that can be obtained in practice. Bet- 
ter results can be obtained from numerical optimization. 
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Rates 


Rates (current) 


Q (current) 


Rates (better) 


Q (better) 


WV\ Direct coupling 


9 


< 100 


< 272 


< 200 


< 619 


ijV] Virtual interaction 




~ 10 


~ 16 


~ 20 


~ 31 


gVIIBI Red /Blue sideband 


g(n Rl /2A ai )(n R2 /2A a2 ) 


> 1 


> 3 


> 2 


> 6 


3 VII CI Red /Blue sideband (ac-dither) 


g(Q R /2A a )J 1 (8E c n a g c /E J ) 


> 3 


> 8 


> 6 


> 19 


ijVIIII Geometric ohase 


> g 2 /200A r 


> 3 


> 4 


> 6 


> 9 



TABLE I: Rates and quality factors of the various gates. For the direct coupling gate, the values presented are upper bounds. 
For the sidebands and geometric phase gates the results presented are only lower bounds. The system parameters are taken as 
g/A = 0.1 = 0k/2A o = 0.1 for the approximations used to derive the rates to be valid. For the 'current' results, we have taken 
g/2-K = 100 MHz and n/2ir = 0.1 MHz. For the 'better' resultSjWe have taken g/2ir = 200 MHz and k/2tt = 0.01 MHz. The 
dephasing time was taken as T2 = 500 ns, as measured in Ref. [3| in this system. For the geometric phase, we took e = A r to 
satisfy both Eq. (|8.18|l and the dispersive approximation. For the "current" results, we have taken A r /2n — 32 MHz and for 
the "better" results A r /2vr = 70 MHz. 



This was already shown in Fig.[7]for the red and the blue 
sideband and also in section lVIIII for the geometric phase. 
Therefore, there is good hope that all of the gates dis- 
cussed in this paper could be realized experimentally, ad- 
mittedly with different degree of usefulness for quantum 
information processing. Moreover, all the quality factors 
quoted in table U arc limited by the system's decoherence 
time T 2 and not by measurement-induced dephasing rate 
from the application of the gate. As a result, increasing 
T-2 in this system will lead to significantly better quality 
factors. 

The key to improve the gate's quality factor is to im- 
prove the coupling coupling g, the resonator and qubit 
relaxation and de pha sing times. Since the first circuit 
QED experiment [15| . g as been improved by almost a 
factor of 20 and further improvements can be realized 
without technical challenges. Resonators with long pho- 
ton lifetimes, in the tens of microseconds, have already 
been fabricated [5(| and first steps in the design and re- 
alization of a new charge-type qubit which is largely in- 
sensitive to charge noise have been taken [l7|. Circuit 
QED therefore seem like a promising system with which 
to study quantum mechanics at the large scales and quan- 
tum information processing. 
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APPENDIX A: FIELD DISPLACEMENT ON THE 
MASTER EQUATION 

In this appendix, we apply the field displacement pro- 
cedure of section IIIII on the the master equation (|2.5p 
rather than on the pure state evolution. Since the dis- 
placement is on the resonator field, we only consider the 
contribution to damping due to k. The master equation 
thus reads 

p = — i[H, p] + — (2apa) — a)ap — pa) a) , (Al) 
where 

H = u> r a) a H — jr* 7 * — 9 (a) <j- + cr + a) 
+ J2 (^(Ve^ 1 + e* k {t)ae +i ^ t ) . 

k 

For simplicity of notation, in this appendix we only con- 
sider a single qubit and drop the j index. Going to the 
displaced frame, the master equation for the displaced 
density matrix p = D' {a)pD{a) reads 

p = - i[D^(a)HD(a),p} - D f (a)D(a)p - pD f (a)D(a) 
+ -£)l(a) (2apa f - a)ap - pa f a) D(a) 

= — i[H, p] + — (2apa) — a) ap — pa) a) 

(A3) 

where 

H = to r a) a H -a r — q (a) O- + o~+a) 

2 yV + ; (A4) 

— g (a*(T_ + ao~ + ) 
and the parameter a is chosen to satisfy 

a = -iuj r a - i'^2eke~ luJd * t - -a. (A5) 

k 

As an example, we consider the simple case of the bit- 
flip gate discussed in section IIII Al In this situation, a 
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single timc-indcpcndcnt and real drive e is needed. With 
a(0) = and dropping the transient term, we recover the 
Hamiltonian (13.51) where the Rabi rate now reads 



R = 2- 



in/2' 



(A6) 



As it should, this rate does not diverge at A r = 0. Since 
we always work in the regime where A r is large, we ne- 
glect this k correction in most sections of this paper. 



APPENDIX B: DERIVATION OF THE 
EFFECTIVE HAMILTONIAN 

In this appendix we derive an effective Hamiltonian for 
two qubits dispcrsively coupled to a single resonator and 
taking into account the presence of three independent 
drives. Two of these drives (k = 1,2) are used to co- 
herently control the qubits while the third drive (k = 3) 
plays the role of measurement beam. Our starting point 
is Hamiltonian (|2.3p taking into account the gate charge 
dependence. While working away from charge degener- 
acy is not useful in practice, this will make selection rules 
appear clearly. 

The derivation follows the same steps as those pre- 
sented in section IIIII For simplicity, we take the drive 
amplitudes Ek to be time- independent. This assumption 
is relaxed in the context of the geometric phase gate in 
section IVIII1 As in section IIIII we start by displacing 
the resonator field a by using the displacement operator 
D(a). Since £3 plays the role of measurement, it's fre- 
quency will be close to the resonator frequency u/ r . It is 
therefore more convenient to displace the field only with 
respect to the first two drives. In the lab frame, the result 
of this transformation is 

ff« = uirJa + -fo Zl + -^a Z2 + e 3 (aV^ + h.c.) 
- X! g J ~ c i a *i + s i ax i) ( " + a ) 

3=1,2 

+ ^R Jk {c 3 a Zj ~s 3 a Xj )cos(ui dk t), 

j,k=l,2 

(Bl) 

where &,R jk = 2g 3 £k/{uj r — u> dk ) is the Rabi frequency of 
qubit j with respect to drive k. 

In the next step, we make the rotating wave approxi- 
mation on the drive terms acting on the qubits [last term 
of Eq. (HQ ]: 



j\fc=l,2 



(B2) 



Not doing this approximation would only lead to small 
Bloch-Siegert shifts in the qubit transition frequency [72j . 



In this appendix, we are not interested in direct transi- 
tions of the qubits (i.e single-qubit Rabi flopping), these 
transitions are therefore eliminated using the transfor- 
mation of Eq. (|3.10p on each qubit: 



Uj = exp [/3*cr +j - /3j(7_J 



(B3) 



Since [Ui, U2] = 0, these two transformations can be ap- 
plied sequentially: 

ff (2) = (U^H^iU^Y - i Ufir 1 . (B4) 

.3=1,2 

We expand the first term to second order in the small 
parameter (3 3 using the Hausdorff formula Eq. (|3 . 8[) . For 
the second term of Eq. (|B4[) . we obtain to second order: 



I / \ ' « ~ m)o Z] + fa*-, - $*a +j ). (B5) 



Choosing 



iuj d t 



(B6) 



fc=l,2 J fc 



with A ajk = u aj — u)d k and neglecting fast oscillating 
terms, we finally obtain 



3=1,2 

~ X 9j {Vj + SjVxj) (a f +a) 
3=1,2 

- X 9jA 3 ((3*<r +j + &cr_.) (a* + a) 
3=1,2 

~ X 9]B 3 a Z] (a) + a) , 
3=1,2 

where the shifted qubit transition frequency is 



(B7) 



= ""J + X 



k=l,2 ajk 

2A ajl A aj2 



(B8) 



cos(w dl - u d2 )t 



and 



A 3 =2c 3 -s 3 {p*+p 3 ) 

B 3 =2c 3 (iPtf-l) +s 3 ((3*+P 3 ). 



(B9) 



In the presence of a single off-resonant drive, we recover 

Eq. (HUD. 

Finally, the qubits are assumed to be strongly detuned 
from the resonator. We therefore adiabatically elimi- 
nate the direct Jaynes-Cummings qubit-rcsonator inter- 
action. This is done using the dispersive transformation 
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of Eq. (|5.1[) . Since the rotating wave approximation was 
not performed on the qubit-rcsonator interaction, this 
choice of transformation will not cancel completely the 
interaction. Complete cancelation could be obtained by 
choosing a slightly different transformation. However, for 
the transitions of interest here, the remaining terms will 
be oscillating rapidly and will simply be dropped. Taking 



into account these terms would, again, only add a small 
frequency shift [72j that can safely be ignored here. 

Applying the dispersive transformation Eq. (|5.ip . ex- 
panding to second order in the small parameter ffj/A^- 
and neglecting fast oscillating terms, we obtain the main 
result of this appendix 



« u, rfl t a + J2 ^fcx Zj + e 3 (ate-*--.* + ae +i ^) - £ f fel (aV^ + ae +^ 3 ') Gz . 

j=l,2 3=1,2 V 3 ) 

sr r M, u giff2(A , 1 + Ag) 

- 2^ 9i B i a *A a +a)+sis 2 K^- 2 + a_ l( 7+ 2 ) 

3=1,2 1 2 

£ '•,'/, K^' + .-,e+^V) ( a t + „) 

j,fe=l,2 ^ 

+ > s j9j~z~ ^ — il e fc' ; +e 1 fc Jer +j . + le v k >" 1 + e v "* V J cr-^. > (a 1 + a), 



where the shifted qubit transition frequency is 

<=<+2s2j-(ya+£) (Bll) 
and we have defined 

A;- = J a . - oj r . (B12) 

This effective Hamiltonian contains all the physics 
needed to realize each of the different gates that are stud- 
ied in this paper. More particularly, the first term in the 
second line of (|B10|) can be used to generate a geomet- 
ric two-qubit phase gate [H, |4?J and is studied in more 
details in section IVIIII The second term of the second 
line of equation (|B10[) is the flip-flop interaction due to 
virtual interaction with the resonator and is discussed in 
section [V] We note that higher order flip-flop terms in- 
duced by the external drives have been dropped. Finally, 
as discussed in section IVH B[ the last two lines describe 
one and two photon blue and red sideband transitions. 

Higher order terms in the perturbative expansions used 
to obtain Eq. (|B10|) will lead to additional non-linear 
terms. These corrections will be negligible as long as the 
parameters f3j and gj/Aj are chosen to be small. This 
can be safely done in the context of the quantum gates 
studied here. The situation where these terms can no 
longer be neglected will be discussed elsewhere. 



APPENDIX C: RED AND BLUE SIDEBANDS 
USING FLICFORQ 

Here we consider a single qubit coupled to the res- 
onator and in the presence of a single drive of frequency 
u) a and amplitude s. For simplicity, we work at charge 
degeneracy and in the rotating wave-approximation: 

H = uj r a) a H -a z — q(a) <J- + u+a) 

2 (CI) 

+ e(^e~ iUat + «?*"•*). 

We again assume that the qubit is far detuned from the 
resonator frequency. However, we choose not to adiabat- 
ically eliminate the qubit-rcsonator interaction and look 
at a non-perturbative result for large drive amplitude. 
To find the effective qubit-resonator Hamiltonian in this 
case, we move to a frame where a non- vanishing interac- 
tion remains, but only when a resonance condition to be 
determined is satisfied. 

We first move to a frame rotating at the qubit transi- 
tion frequency oj a for both the qubit and the resonator: 

H (1) = -Aa'a - ff(aV_ + a + a) + e{a) + a), (C2) 

with A = uj a — oj r . We then displace the resonator field 
using the displacement operator D(a) to obtain 

H {2) = -Aft+fl - g(a}a- + a+a) - ^y<?x- (C3) 

For convenience, we now change a x to — a z using a rota- 
tion along the y axis of angle —ir/2. After that rotation, 
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the system looks like a qubit of frequency £Ir coupled to 
a resonator of frequency A: 



= -Aa+a + -^a z + |(o+ + a)a z + i|(o+ - o)a tf . 

(C4) 

Going to a frame rotating at both the new effective qubit 
and cavity frequencies, we obtain a single non-oscillating 
term when A = —Qr: 



(aV 



and for A = £Ir we obtain: 



(4) _ 



t 

— ((2 (7^. 



.a). 



.a). 



(C5) 



(C6) 



These correspond, in the rotating frame, to red and blue 
sideband transitions at a rate g/2. The factor of 1/2 is 
due to the fact that only one of the Rabi sidebands for the 
ground and excited states of the qubit are in resonance 
at A = This type of qubit-cavity resonance was 

also discussed in Refs. 
a dressed basis in Rcf 



73 1 and 7J| and for flux qubits in 
19 . 



APPENDIX D: FLICFORQ WITH TWO QUBITS 

To obtain the effective Hamiltonian for the two-qubit 
FLICFORQ, we follow the results of appendix iBl and of 
Ref. [5^]. The starting point is the Hamiltonian (jBip 
taken in the rotating wave approximation and at charge 



degeneracy for both qubits (sj 



l,Ci 



0). Moreover, 



we omit the measurement drive and consider only two co- 
herent drives (k = 1, 2). To get maximal splitting of the 
Rabi sidebands, the frequency of these drive are chosen 
such that LOdj is close to cu aj . 

To derive the effective Hamiltonian, we follow the same 
steps as in appcndix|B] The main difference is that in the 
second step, because of our choice of drive frequencies, 
we do not adiabatically eliminate both drives from both 
qubits. We adiabatically eliminate only the effect of co^ 
on the second qubit and the effect of Ud 2 on the first 
qubit. The resulting effective Hamitonian is 



H m uo r a) a + — -cr Zj 



3=1,2 



+ 



3=1,2 
5152 (A' x + A2) 



(Dl) 



2A;A' 2 



where A' 



oj r and 



Ol(2) 



, -"12(21) 
^ a 12(21) 



(D2) 



Assuming for simplicity that — , we first go to 
a frame rotating at ojdj for qubit j. Following ref. [59l ]. 



we then go to a frame rotating along the x axis and at a 
frequency Q,rJ2 for qubit j. This yields 



tt~ t , 9m{A[ +A' 2 ) 
ti ss uj r a' a + 



8AiA^ 

x {Wxt - i sxn(Cl Rl t)a Zl + icos{£l Rl t)a yi ] 



x cr :r 



i sin(0fl 2 t)cr Z2 — icos(Q,R 2 t)a. 



!J2\ 



-h.c.}. 



(D3) 



We now look for terms that do not average to zero. These 
will correspond to resonance in the coupled driven sys- 
tem. There are several choices of resonances here and, as 
an example, we choose u>d 2 — ^>di = fifli + Qr 2 - Using 
this condition, all terms but a single term are oscillating 
rapidly and average out to zero. The resulting effective 
Hamiltonian in the quadruply-rotating frame is 



TT t , 3152 (Ai 

-Hff ~ uj r a' a+ 



16AiA^ 



,)■ (D4) 



Other choices of resonances lead to different symmetry 
for the effective Hamiltonian. These other possible cou- 
pling Hamiltonians will be discussed elsewhere. 



APPENDIX E: SYMMETRY AND SELECTION 
RULES 



At the charge degeneracy point, the Hamiltonian (12.31) 
is even in the number of creation and annihilation oper- 
ators. As a result, C = a^a + a z /2 commutes with the 
Hamiltonian and the total excitation number is a con- 
served quantity. On the other hand, away from charge 
degeneracy, in addition to the regular Jaynes-Cummings 
term, we get extra terms of the form g(fi + cer 2 )(a^ + a) 
which arc of odd parity in the number of creation and an- 
nihilation operators. Clearly, C is not a conserved quan- 
tity in this case. 

These symmetry considerations will impose selection 
rules on the transitions that are possible. To see this 
more explicitly, we introduce the parity operator P = 
e~ ma a a z which is the natural unitary extension of C [75l . 
[76j . By writing the parity operator in the form P = 
X)^Ln( — l)"l n ) ( n Wz, h is simple to verify that P anti- 
commutes with the drive Hamiltonian (|3.ip in a frame 
rotating at the drive frequency. Labeling states as being 
of even (odd) parity if they are eigenstates of P with +1 
(-1) eigenvalue, it is then clear that the drive (al + a) can 
only cause transitions between states of different parities. 

To see which transitions are allowed at the sweet spot, 
we first give the parity of the eigenstates of the Jaynes- 
Cummings Hamiltonian: 



P\g0) 

where [lj] 



■\g0), 



P\±^) = (-l) n \^n), (El) 



6>„|e,n) +sin6 l „|g,n+ 1) 



|+,n) = cos 

| — ,n) = — sin6*„|e, n) + cos n \g, n 



(E2) 
(E3) 
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and 



1 , -if 2gvn 
— tan 
2 



(E4) 



The red sideband transition illustrated in Fig. [6^,) cor- 
responds to a transition between |+,0) and | — ,0), it is 
therefore forbidden to first order. This is the result al- 
ready obtained in Eq. (|7.5p . In the same way, the blue 
sideband corresponds to a transition between \g0) and 
| — , 1) . It is also forbidden at charge degeneracy. Single- 
photon Rabi flopping discussed in section fill Al connects 



\g0) and | — ,0), and it is obviously allowed at charge de- 
generacy. In the case of two photon transitions, the drive 
Hamiltonian is effectively acting twice and the previous 
selection rules are therefore simply reversed. Finally, 
away from the charge degeneracy point, the Hamiltonian 
has no definite parity and these considerations do not 
apply. All sideband transitions are therefore allowed. 
Similar selection rules for flu x q ubits cou pled to a LC 
oscillator were studied in Ref. [2fJ| . In Ref. [70| , these se- 
lection rules were investigated for a flux qubit irradiated 
with classical microwave signal. 
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